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■ ABSTRACT 

We study the properties of Nambu monopoles and Z-vortices in 
^ . the ?>D lattice SU{2) Higgs theory which represents the Standard 

Model at high temperature. We show that the densities of the Nambu 
\ monopoles and the Z- vortices are 0(1) in the symmetric phase and 

^ \ generically small in the Higgs phase. Near to the critical Higgs mass 

and in the vicinity of the phase transition the densities are no more 
negligible in the broken phase. The percolation probability of the 
Z-vortex lines is found as a new disorder parameter for this phase 
transition. We conclude that the transition to the symmetric phase is 
accompanied by Z-vortex condensation. Simulations comparing ele- 
mentary and extended vortices and monopoles at different /3g values, 
aiming to show that the density of vortices and monopoles of fixed 
physical size might have a well-defined continuum limit, gives encour- 
aging but so far inconclusive results. 
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1 Introduction 



This paper is devoted to the study of topological properties of thermal states in the 
standard model of electroweak interactions. Although the standard model does not possess 
topologically stable monopole- and vortex-like structures, one can define so-called embedded 
topological defects [|I], : they are known under the names of Nambu monopoles and 
Z-strings 0, We perform the investigation of the equilibrium thermodynamics of 
embedded topological defects in the framework of dimensional reduction which is expected 
to be reliable to describe the physics near to the transition temperatures (of a truly first 
order phase transition or a continuous crossover) for Higgs masses between 30 and 240 
GeV. 

Embedded topological defects play a crucial role in some electroweak baryogenesis 
scenarios. For example, one of the scenarios is based on the decay of an electroweak 
string network as the Universe cools down. According to this mechanism, long elec- 
troweak strings decay into smaller, twisted and linked string loops which carry non-zero 
baryon number. This could then explain the emergence of non-zero baryon number in 
the Universe |^. Another example is the bound state picture of a sphaleron: a Nambu 
monopole-anti-monopole pair connected by a segment of Z-string p, |^, ^ builds the 
unstable (saddle point) sphaleron field configuration ^ in the broken phase. After 
the electroweak transition is completed, this type of electroweak configuration is the only 
transition state available between vacua of different Chern-Simons number. Thermal ac- 
tivation would allow to wash-out the baryon number asymmetry and therefore has 
to be suppressed. In a recent paper we have analyzed classical lattice sphalerons (as 
obtained and discussed in Ref. |T0[) and have found evidence for the monopole-anti- 
monopole bound state picture of these configurations. 

In the present publication, we apply the tools for detecting Nambu monopoles and 
Z-vortices developed in Ref. ^ in order to study numerically the formation of these 
defects in thermal equilibrium states of the electroweak theory. We measure the densities 
of Nambu monopoles and of Z-strings across the phase transition. In an attempt to sub- 
stantiate the picture of Z-vortex condensation we also study the percolation probability 
of the Z-vortex line network near to the phase transition temperature. Recently, in a dif- 
ferent context, in Ref. [|TT| convincing evidence has been provided that string percolation 



is a good disorder parameter for a phase transition in field theory. 

Due to the smallness of the Weinberg angle 6w the U{1) component of the electroweak 
group SU{2) X f/(l) can have only little effect on the properties of the embedded defects. 
Thus the dynamics of the topological defects can be studied restricting ourselves to the 
SU (2) Higgs sector. The emerging effective 3D SU{2) Higgs lattice modelQ allows to 
take into account effects of chiral fermions on physical parameters perturbatively [|13|, |14 . 



"Besides symmetric and Higgs phases, the lattice SU (2) Higgs model in AD has also the confinement 
phase which is located at sufficiently small Pg (which has, however, nothing to do with the electroweak 
theory in the continuum). One can show analytically that the transition between the symmetric and the 
confinement phase in that model is accompanied by Nambu monopole condensation . This observation 
leads us to conjecture that topological defects may be relevant for the phase transitions also in the 
physically relevant region of the parameters of the SU{2) Higgs model, which is studied here in its 
dimensionally reduced variant. 
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However, as in our study of the structure of the electroweak sphaleron one may expect 
that the influence of the standard model fermions on the embedded defects is small [Q. 

The structure of the paper is as follows. The effective three-dimensional model is 
shortly recalled in Section 2. In this Section we also give the lattice definitions of the 
embedded topological defects according to Ref. and generalize them to cover the case 
of extended defects, too. In Section 3 we present the numerical results on the density 
of elementary (size a) monopoles and vortices and on the percolation probability of the 
corresponding Z-vortex lines. In Section 4 we describe a first step to extend this analysis 
to extended vortices and monopoles of equal physical size living on finer lattices in order to 
see whether the density discontinuities at the transition temperature possess a continuum 
limit. Section 5 contains a discussion of our results and conclusions. 



2 Embedded Defects in Lattice SU{2) Higgs Model 

We study the lattice 3D SU (2) Higgs model with the following action: 

S = /3c5^(l-lTr[/,)-/?^^iTr(<|.+f/.,,$.+^) + 5:(p^. + /5^^^^^ (1) 

where the summation is taken over plaquettes p, sites x and links / = {x, /i}. The action 
contains three parameters: the gauge coupling /^g, the lattice Higgs self-coupling and 
the hopping parameter ^h- The gauge fields are represented by unitary 2x2 link matrices 
Ux,ii. and Up denotes the SU{2) plaquette matrix. The Higgs fields are parametrized 
as follows: = PxVx, where pi = |Tr($+<l>2,) is the Higgs modulus squared, and 
an element of the group SU{2). The lattice parameters are related to the continuum 
parameters of the 3D superrenormalizable SU{2) Higgs model g^, A3 and ms^fj,^ = g1) as 
given e.g. in |[T^. As in ||14| a parameter is used (approximately equal to the zero 
temperature physical Higgs mass) which expresses the Higgs self-coupling as follows 

gUc sUoGeV^ ^g' ^ ' 

The lattice coupling Pg and continuum coupling (7I are related by 

Pg = ^, (3) 

with a being the lattice spacing. We shall study the theory at different gauge couplings Pg 
in order to understand the qualitative behavior of the lattice embedded defects towards 
the continuum limit, i.e. when the lattice becomes fine enough to resolve the structure 
of defects which then have to be described as extended defects. 

Let us first define elementary topological defects. The gauge invariant and quantized 
lattice definition of the Nambu monopole is closely related to the definition in the con- 
tinuum theory ||^. First we define a composite adjoint unit vector field = n°^(T°', 

= — ($+o""$^.)/($+$^) with a"" being the Pauli matrix. The field n^. plays a role sim- 
ilar to the direction of the adjoint Higgs field in the definition of the 't Hooft-Polyakov 
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monopole |T5| in the Georgi-Glashow model. Then we define the gauge invariant fiux 
plaquette 9p: 

e,{U,n) = arg(Tr [(1 + n.)\4,^\4+A,.V;+ ,,^V;+ ]) . (4) 

Here 

Vx,^,{U, n) = U^^^ + nJJ^^^,n^+f, (5) 

has been used to define a particular projection of the gauge field. The sense of it becomes 
more transparent if we consider the unitary gauge where we have $2; = (0, 0)"^, = 5°"^ . 
In this gauge the phases 6*" = arg f//^ behave as a compact Abelian field with respect to 
the residual Abelian gauge transformations VL°^^^ = e*'^^""', ax € [0,2-7?) which leave the 
unitary gauge condition intact. The Nambu monopoles are the topological defects in this 
Abelian field. The gauge invariant plaquette function coincides with the standard 
Abelian plaquette formed out of the fields 6f [0]. In this case, the Nambu monopole 
charges can be defined using the standard construction: 

jc = -^Y^Op, ep = {Bp - 2TTmp) G [-TT, tt) . (6) 

In the gauge independent definition (Q) the fiux immediately gets values in the required 
interval [— tt, tt) and is independent of the choice of the reference point x on the plaquette, 
due to the property n^K,/. = K./^^-x+a- 

The Z-string |^ ^ corresponds to the Abrikosov-Nielsen-Olesen vortex solution [|r^ 
embedded Q into the electroweak theory. The Z-vorticity number corresponding to 
the plaquette p = {x, /iz/} is defined similar to the construction of the Abelian vortices 
in an Abelian Higgs theory: 

= ^{Xp - dp) , (7) 

where 6p has been given in eq. (D, and Xx,,iu = Xx,t, + Xx+ti,u-Xx+u,fi-Xx,u is the plaquette 
variable formed out of Abelian links 



Xx,^^ = arg ($+K,A.<^'x+/i 

The Z-vortex is formed by the links / = {x, p} of the dual lattice carrying non- 
vanishing vorticity *ai. These links are dual to the plaquettes p = {x,iii>} with non-zero 
vortex number (|^ and *<Jx,p = £p^lv0^x,^lv/'2,. One can show that Z-vortices begin and end 
on the Nambu (anti-) monopoles: E^=i(*o-x-/i,/. - *CFx,p) = *jx- 

In order to understand the behavior of the embedded defects towards the continuum 
limit we have made measurements of densities of Nambu monopoles and Z-vortices just 
of one definite physical size which amounts to different multiples of a on different lattices. 
We follow the change from the broken to the symmetric phase at various gauge couplings, 
Pg = 8, 16,24, i.e. with a lattice step size a getting smaller with increasing Pg- In prin- 
ciple, starting from a lattice simulated at various jSc, we could perform different factor k 
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real space renormalization group transformations of the gauge-Higgs field configurations 
in order to measure properties of defects of many different sizes ka objects in a manner 
described above. Instead of doing this, as a first step, we use here the notion of extended 



topological objects on the lattice as invented in Ref. [17|. An extended monopole (vortex) 
of physical size ka is defined on k^ cubes (fc^ plaquettes, respectively) .Q The charge of 
monopoles jc(fc) on bigger k^ cubes c{k) is constructed analogously to that of the elemen- 
tary monopole, eq. (^, with the elementary 1x1 plaquettes in terms of V^^^ replaced by 
n X n Wilson loops (extended plaquettes) denoted as p{n). 

First, the projected links Vx^^{U,n) are defined as usual, then extended \^-links 

^xji = yx,tJ.Vx+fi,f,...Vx+{k-i)fi,,i (9) 



are formed where X and X + kfi denote points of the decimated lattice, which inherit 

(k) (k) 

the property nxVx/f^ = Vx/f^nx+kji- Finally, the phase angles of length-fc extended Higgs 
operators are defined as 

Xm = X^Z = arg • (10) 

With the gauge invariant flux ^p(fc) for the extended plaquettes defined analogously to 
(I) in terms of the extended y-links and with the extended plaquette angles Xp{k) formed 
out of the extended links angles xz(fc) , the definitions of extended monopoles and vortices 
can be taken over from eqs. (|],^ 

jc{k) = -TTI H ^p(fc) > f^p(fc) = ^{Xp{k) - ^p(fc)) • (11) 

p{k)€dc{k) '^^ 

According to eq. (§) the physical size of the k^ monopoles (or the fc^ vortices) in simu- 
lations done at /?g = ^/^g"* should be equal for all k. Roughly the same physical objects 
are hopefully described as elementary monopoles and vortices at i3q\ too. 



3 Scanning the Transition with Size a Defects 

First, the formation of elementary Nambu monopoles and vortices has been scanned 
across the electroweak phase transition. Monte Carlo simulations have been performed 
on cubic lattices of size = 16^ at Pc = 12 for self-couplings A3 corresponding to 
M|f = 30 GeV and = 70 GeV, see eq. (H). We varied the parameter Ph in order to 
locate the phase transition for given and (3g- For the update we used the algorithms 

''In the context of Abelian projections of pure Yang-Mills theory this construction is known as type-I 
extended monopoles So-called type-II monopoles have been defined in the same context. These are 
defined by adding the charges (currents) of elementary monopoles inside the extended cube c(n), and 
type-II vortices can be constructed similarly. We have also tried this construction for our model and 
found that here this block spin transformation does not seem to lead to a well-defined continuum limit 
for extended defects contrary to the case of SU{2) gluodynamics where the type-II Abelian monopoles 
behave according to the renormalization group [ psf . In the limit a — > the detection of size a topological 
defects with the definition given above is afflicted with increasing lattice artifacts. 
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described in Ref. which combine Gaussian heat bath updates for the gauge and Higgs 
fields with several refiections for the fields to reduce the autocorrelations. Usually, the 
first order nature of the thermal transition is detected by a two-peak histogram in the 
volume average of p^. The average over all Monte Carlo configurations (p^) denotes the 
quadratic Higgs condensate, the discontinuity (for infinite volume) of which is related to 



the latent heat 14 



Here we are interested in the behavior of the lattice densities pm of Nambu monopoles 
and Pi, of Z-vortices and want to measure the percolation probability for the Z-vortex 
lines near to the phase transition. For each lattice configuration, the densities pm and p^ 
are given by 

P^ = -kY.\3c\. P^ = ^El^pl' (12) 

where c and p refers to elementary cubes and plaquettes. The monopole charge jc and 
the Z-vorticity ap are defined in eq. (|) and eq. (|^, respectively. 

The percolation probability of the vortex lines *a is defined with the help of the 
following two-point function 



C(^) = II<^xe*aW'^ye-aW-'^(k-Z/|-^) ■ I]5(|a;-l/|-r) , (13) 




where the summation is taken over all points x, y of the dual lattice and over all connected 
clusters of vortex lines *cr*-*^ {i labels distinct vortex clusters). The Euclidean distance 
between two points x and y is denoted as \x — y\. The notation x G *a^'^'^ means that the 
vortex world line cluster *a'^'^^ passes through the point x. 
Then the percolation probability is 

C = lim C(r). (14) 

r— ►oo 

This formula corresponds to the thermodynamical limit. In this case the function C(r) 
can be fitted as C(r) = C + Cor~°'e~''^^ , with C, Co, a and m being fitting parameters. 
In our finite volume we find numerically that m ~ in the studied region of the phase 
diagram. Therefore we can be sure that finite size corrections to C are exponentially 
suppressed. If C does not vanish, i.e. if two infinitely separated points with non-zero 
probability belong to the same connected cluster of vortex lines, the vacuum is populated 
by infinitely long vortex lines. In other words, an effective theory of these vortices would 
have an action smaller than the entropy per length of the vortex trajectory. This case 
implies the existence of a vortex condensate. If C turns to zero the vortex condensate 



vanishes 19 



In Figure 1(a) we show the ensemble averages of densities, {pm) of Nambu monopoles 
and (p^) of Z-vortices vs. hopping parameter Ph for Higgs mass = 30 GeV at 
gauge coupling (3g = 12. At this Higgs mass the phase transition is of strongly first 
order [^. It is seen that the averages of both densities are 0(1) in the symmetric phase 
and they sharply drop at the phase transition to zero densities in the Higgs phase. In 
Figure 1(b) the percolation probability C for the same parameters is shown to have a 
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small discontinuity at the transition. Deeper into the symmetric phase the condensate 
(percolation probability) grows while it is strictly vanishing in the broken phase. 

For comparison, Figure 2(a) presents the Monte Carlo averages (p^) and (p^) for the 
Higgs mass = 70 GeV at jSc = 12. This Higgs mass parameter is very near to the 
critical Higgs mass. Now, the average densities of Nambu monopoles and Z-strings do not 
vanish immediately after the phase transition towards the broken phase: both densities are 
non-zero at temperatures T <, T^. Still, the corresponding vortices do not percolate. The 
percolation probability C for the Z-vortex trajectories behaves as a disorder parameter 
also in this case when the high-temperature electroweak phase transition is already very 



weakly first order |T3|. It turns to zero exactly at (3hc as can be seen in Figure 2(b) but 
a continuous behavior cannot be excluded from these data. 

At the larger Higgs mass = 70 GeV, just for the value of the gauge coupling 
(3g = 12 studied before, we have also tried to understand how the infinite volume limit 
is approached. We show in Fig. 3(a) and Fig. 3(b) histograms of the Nambu monopole 
density pm and of the Z-vortices p^, respectively, for lattice sizes 16'^, 24^, 32^ and 40'^. 
Even near to the end of the first order phase transition, we see that not exceedingly 
large lattices are necessary to see the two-state signal in the densities and to measure 
the (finite volume) discontinuity of the densities, such that its infinite volume limit can 
be obtained. In the present state of our analysis we cannot yet present corresponding 
histograms of the percolation probability which would allow to tell whether the weakly 
first order phase transition is still accompanied by a discontinuous change of the vortex 
condensate. Independent of this, one may say that an infinite percolating network is 
breaking into finite pieces at the transition temperature while the correlation length stays 
finite. 



4 Is There a Continuum Limit ? 

In Ref. lattice vortices have been studied in the non-compact version of the 3D Abelian 
Higgs model. The 3D Abelian Higgs model has become a model system for many phase 
transitions, covering the electroweak and the QCD deconfinement transition, supercon- 
ductivity as well as cosmic strings. From their attempt to define a well-defined continuum 
limit of the vortex density two lessons can be drawn. First, the lattice density might be 
afflicted by regularization artifacts such that, generally speaking, discontinuities at the 
phase transition are physically better defined. Second, one should carefully associate the 
density of defects with (various) physical scales and try to find the continuum limit of 
these. At least at the strongly first order phase transition at small Higgs mass the first 
problem seems to be not to so severe. But in order to get a physically well-defined vortex 
and monopole density, the extrapolation to the continuum limit requires the measurement 
of more and more extended defects on finer and finer lattices. 

For a preliminary exploration we have concentrated on the case of Higgs mass = 30 
GeV where the phase transition is of strongly first order. At this Higgs mass the finite 
volume effects are numerically small and therefore we felt safe to study only one lattice 
size, 16^, at various Pg- For the coarsest lattice we have chosen i3g = (3g^ = 8. Just to be 
specific, we have chosen the lattice spacing corresponding to this gauge coupling as the 
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size of topological objects to be examined for a reasonable continuum limit. Now we have 
measured at (3g — k^^Q (for /c = 2, 3) histograms for the densities of extended monopoles 
and vortices of a size (in lattice units) of and A;^, respectively. The measurements were 
done near to the pseudocritical (5h for the appropriate gauge coupling at this lattice size. 

The distributions of Nambu monopole and vortex densities (for objects of corre- 
sponding physical size) are shown in Fig. 4(a) and (b), respectively. These (normalized) 
distributions of extended defects have the same two-peak structure as elementary ones. 
With decreasing inverse gauge coupling 1/(5g the metastability becomes weaker if the 
density of the defects is measured. At the smallest jSc = 8, for instance, where we looked 
for size a defects, we could not observe tunnelling between the two phases within 10000 
iterations following hot/cold starts. In the Figures, for that the histograms corre- 
sponding to the Higgs phase are omitted since only the bin at zero density is populated. 
The peaks corresponding to the symmetric phase shift considerably from elementary ob- 
jects at /Sg = 8 to size-2 objects at (5g — 16, but only slightly from this measurement to 
size-3 objects at Pg = 24. 

We have also checked the percolation of extended vortices. At /3g = 16 on a 16^ lattice 
we have found I3hc — 0.3391. Considering all possible k — 2 decimations, we have now 
analyzed networks of A; = 2 vortices on this lattice with respect to percolation. There 
arc 2^ different k = 2 vortex configurations per gauge-Higgs configuration which have 
been obtained at (3h = 0.3390 and = 0.3392. At Ph = 0.3390 we have found a non- 
zero condensate also for these extended objects, whereas at /3h — 0.3392 the condensate 
is strictly vanishing. We have to admit that this result has been obtained on a very 
small lattice. Nevertheless it is indicative for extended vortices percolating at the same 
temperature. 

Similar measurements have to be repeated for defects of different scale on lattices 
with different lattice spacings and with bigger, physically comparable volumes before 
an understanding of the topological mechanism can be achieved. Notice that the above 
results were obtained concentrating, somewhat arbitrarily, on topological objects with 
smallest possible size that could be analysed on the coarsest, /3g = 8 lattice. Obviously, 
the next task is to study non-elementary topological objects already on the coarsest 
lattice. This endeavour requires much more computer power in view of the complementing 
measurements that have to be done on finer lattices covering the same volume. Since our 
studies already indicate that (within our lattice construction) Nambu monopoles and 
Z-vortices as topological excitations of definite physical size become more and more well- 
defined in the continuum limit, it seems to be preferable to dispense with the elementary 
objects completely. 

The observed, more well-defined behavior is related to the way how the fundamental 
Higgs field is turned into an auxiliary composite field used to project out the Abelian 
fluxes. It seems, that their densities possess a sensible continuum limit which can be 
measured by our lattice prescription. In contrast to this, elementary monopoles and 
vortices are still far from the continuum limit for defects of corresponding size. A similar 
problem is well-known for instantons in 4D lattice gauge theory. 

As indicated above, there is another possibihty to deflne extended defects of any size for 
a sufficiently fine lattice to start with. This requires to perform block spin transformations 
of Higgs and gauge fields which give a consistent description of the phase transition if seen 
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from different length scales. In Ref. |22] two of us have indicated a practical way how to 
do this. 



5 Discussion and Conclusions 

We have studied numerically the properties of Nambu monopoles and Z-vortices which 
appear as embedded defects in the Standard model at high temperatures. Our results 
are obtained in the 3D lattice formulation of the SU{2) Higgs sector. We have shown 
that these defects are dense in the symmetric phase and that they exist as a more or 
less dilute gas in the Higgs phase. Moreover, the percolation probability of the Z-vortex 
trajectories has been shown to be the disorder parameter for the hot electroweak phase. 
The Z-strings are condensed in the symmetric phase, and the condensate vanishes in the 
broken phase. This result has been obtained for the Higgs mass M'^ = 30 GeV at which 
the phase transition is of strongly first order. The same is found to be true for the Higgs 
mass M'^ = 70 GeV (weakly first order phase transition) where the density of monopoles 
and vortices is clearly discriminating the phases. Note that in this case at the phase 
transition the densities in the broken phase are not less than 50 % of the corresponding 
ones in the symmetric phase^. 

Comparing the densities of the embedded defects of corresponding physical size at 
various values of the gauge coupling (3g gives a first hint that the densities of extended 
Nambu monopoles and Z-strings as defined in this paper might possess a well-defined 
continuum limit. Further studies are necessary to understand the dynamics of defects 
and the spatial structure of vortices of different physical size near to the transition and 
in the high temperature phase. 

The approach of this investigation may be extended to lattice SU{2) Higgs model in 
D = 3 + 1. In this case the Z- vortex world trajectories are sweeping out two dimensional 
world surfaces which are open on the monopole world lines. According to analytical 



estimates the Z-vortices form a dense network of long strings at temperatures above 
a temperature which has been termed Hagedorn temperature Tfj in analogy to the dual 
resonance theory of hadronic strings with its exponentially rising mass spectrum. One 
would expect this temperature to be of the order of the critical temperature or slightly 
smaller, <, T^. Below Th the network decays into small loops eventually producing 
non-zero baryon number . We hope to be able to substantiate this picture by AD lattice 
studies in near future. 
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Figure 1: (a) Density of Nambu monopoles pm and of Z-vortices vs. hopping parameter 
13h for Higgs mass = 30 GeV at gauge coupling Pg = 12; (b) Percolation probability 
C of Z-vortex trajectories for the same parameters; the critical hopping parameter is 
Phc ~ 0.3411. 



0.25 



0.20 



>. 0. 1 5 



E 

a 



0.10 



0.05 



o o 



A o 
A 
A 

A° 



A 
b 



Ad 
aO 

I I tii 



^ P. 





0.330 0.335 0.340 0.345 0.350 



0.25 



0.20- 



0.15- 



0.10 



0.05 




0.330 0.335 0.340 0.345 0.350 



(a) 

Figure 2: Same as in Figure 1 for 
Phc ~ 0.34355. 
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Figure 3: Density distributions of Nambu monopoles (a) and of Z-vortices (b) of elemen- 
tary size at pseudocriticality for different lattice sizes, — 70 GeV and Pa — 12. 
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Figure 4: Density distributions of Nambu monopoles (a) and of Z-vortices (b) of fixed 
physical size at pseudocriticality for different gauge couplings, Higgs mass = 30 GeV 
and 16^ lattice. 
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